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I use the two-step density-matrix renormalization group method to extract the critical exponents 
P and v in the transition from a Neel Q — (n, n) phase to a magnetically disordered phase with a 
spin gap. I find that the exponent /3 computed from the magnetic side of the transition is consistent 
with that of the classical Heisenberg model, but not the exponent zv computed from the disordered 
side. I also show the contrast between integer and half-integer spin cases. 



There is current interest in studying quantum phase 
transitions (QPT) in frustrated quantum systems 
An interesting phenomenon generated by frustration is 
the dimensional reduction recently observed in the Bosc- 
Einstein condensation QPT of BaCuSiiO§ [2j]. In ad- 
dition to their connection to real materials, these 2D 
models could be an example of systems where the con- 
ventional Landau-Ginzburg- Wilson (LGW) approach to 
phase transitions may fail. Of particular interest is the 
transition from an antiferromagnet (AFM) to a disor- 
dered phase with a spin gap and short-range magnetic 
correlations known as a valence bond crystal (VBC). A 
very fruitful approach to the study of these transitions 
relies on to the mapping of the Heisenberg model to 
the non-linear sigma (NLct) modeli, Hfl. If the Berry 
phase terms are neglected at the transition as suggested 
in Ref . 0] , the critical behavior of the quantum 2D model 
is identical to that of the NLcr model. However, the Berry 
phase effects are absent only in the AFM phase; they are 
important in the disordered phase and they might well 
modify the critical exponents of the 2D Heisenberg model 
as suggested in Ref.[10|. 

The transition from an AFM to a spin gap in the 2D 
Heisenberg model was studied in Ref. [3J] by the Quantum 
Monte Carlo (QMC) method. It was concluded that this 
transition belongs to the universality class of the 0(3) 
classical Heisenberg model as may be expected from the 
NLa mapping. In this study, the disordered phase was 
generated by the explicit dimerization of the bonds. This 
would support the view of Ref. that Berry phase effects 
do not affect the transition. But due to the minus sign 
problem, the important case of transition driven by frus- 
tration cannot be studied by the QMC method. Hence, 
it remains to be seen whether frustration merely brings 
some technical complication, as would suggest the hy- 
pothesis that the universality class does not depend on 
the detail of the Hamiltonian, or it can drive the transi- 
tion to another universality class. 

There is also a more general interest in studying frus- 
tration induced transitions. The conventional approach 
to quantum phase transitions is to use the path inte- 
gral formalism which maps quantum phase transitions in 
dimension d to classical transitions in dimension d + 1 



to which the machinery of the Landau-Ginzburg- Wilson 
theory is then applied. However, this mapping is fraught 
with difficulty in some important situations. For some 
experimentally relevant quantum mechanical models, the 
corresponding classical functional integrals have non- 
positive or even complex valued Boltzmann weights. This 
would suggest in such cases that the quantum transition 
does not have a classical equivalent. It is thus important 
to study the quantum models directly. But the direct 
study of such quantum mechanical models has proven 
to be a formidable challenge to condensed-matter the- 
orists both numerically and analytically. In particular, 
the QMC method, which in its formulation uses this 
quantum-classical mapping, is at present unable to ad- 
dress these issues because of the non-positive values of 
the effective Boltzmann weights. This state of matter 
has stalled progress in the understanding of QPT. 

In a recent publication [J], it has been shown that 2D 
problems may be attacked by a chain perturbation the- 
ory method. This method called the two-step DMRG 
has recently been applied to study quantum phase tran- 
sitions induced by frustration in two dimensions. The 
two-step DMRG is not a simple Green's function per- 
turbation expansion which is known to fail in situations 
where quantum fluctuations are important. It is rather 
a perturbation expansion on the reduced Hamiltonian. 
It thus retains the full low-energy many-body dynam- 
ics of the original problem. The two-step method has 
the ability to reach an ordered phase when it exists. In 
this letter, we report the computation of the critical ex- 
ponents in the important case of the transition from an 
AFM to a VBC induced by frustration for 5=1. Such a 
transition was suggested in Ref. [13] from a large N anal- 
ysis and found in an anisotropic 2D Heisenberg model 
with 5 = 1 in Ref. 4] . In addition we highlight the differ- 
ence between integer spin and half-integer spin systems 
by comparing the case of 5 = § to that of 5 = 1 . 

I study the following Heisenberg model on the 
anisotropic square lattice: 
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where S = 1, Jii is the in-chain exchange parameter and 
is set to 1; Jj_ and J a are respectively the transverse and 
diagonal interchain exchanges. 

In the TSDMRG, we start by applying the DMRG 
to a single chain for which we obtain the low energy 
eigenvalues e„ and eigenvectors | </>«). Then the Hamil- 
tonian ([1]) is projected unto the tensor product I^iim) = 
l^m)!*^) ■■■\4 , ni) yielding the effective low energy Hamil- 
tonian 

H «^^|| W |$|| W )($|| W | + J ± ^S. M S M+1 + 

[n] il 

Jd Si^Si+i^+i + Si+i^S^z+i, (2) 

a 

where Siir n i is the sum of eigenvalues of the different 



chains, E 



are the renormalized ma- 



trix elements in the single chain basis. They are given 
by S" ! ; '" i! = (<f> ni \Sn\4> mi ). Since the diagonalization has 
been made in the direction of the chains, the effective 
Hamiltonian([2]) is ID. I again use the DMRG to obtain 
its spectrum. 

The TSDMRG has been extensively checkedjl @], and 
is variational. It is controlled by two parameters mi 
and m-2 which are the numbers of states kept during the 
first and second step respectively, mi should be large 
enough so that not only the ground-state energy of a sin- 
gle chain / but all the e ni and <j) ni retained to be used 
during the second step are accurate enough. And, m 2 
should be large enough so that J±, Jd <C e ni ~ e 0i • Fot 
instance for S = 4, we find that the ground-state en- 
ergy per site of a 16 x 17 system evolves from —0.43481 
when (mi, m 2 ) = (128,80) to -0.43681 when (mi,m 2 ) = 
(256, 96) these energies are to be compared to the QMC 
energy —0.43529. It may be argued that despite this ac- 
curacy in the ground-state energy, the TSDMRG would 
not couple the chains effectively and thus may essen- 
tially retain ID physics. This would result in a rapid 
decay of the tranverse spin-spin correlation function C± . 
To counter this argument we show in Table[T] C± for the 
16 x 17 systems. When mi and m 2 are large enough, 
there is no spurious decay of the correlation of the TS- 
DMRG, they appear to decay even slower than the QMC 
ones. Hence, although it starts from an isolated chain, 
the TSDMRG is able to describe the 2D regime. 

In Ref([|), it was found that as Jd moves towards the 
maximally frustrated point Jd ~ Jj_/2 the 2D system 
progressively relax to nearly disconnected chains. At this 
point, the ground state energy is nearly equal to that of 
disconnected chains, the transverse spin-spin correlations 
decay exponentially, and the transverse bond-strength is 
equal to zero up to the numerical accuracy. Hence, for in- 
teger spins, the nature of the ground state is determined 



-0.02116 
0.00726 
-0.00320 
0.00147 
-0.00078 
0.00030 
-0.00013 
0.00006 



-0.03022 
0.01088 
-0.00572 
0.00320 
-0.00186 
0.00108 
-0.00063 
0.00036 



-0.02533(1) 
0.00854(1) 
-0.00399 
0.00201 
-0.00105 
0.00056 
-0.00030 
0.00015 



TABLE I: Transverse correlation C±(l) for a 16 x 17 system 
with J± — 1 and Jd = of two sets of TSDMRG parameters 
against QMC. The origin is taken at the middle of the lattice. 



by the competition between J e f f = J± — 2Jd which favors 
a magnetic phase and the ID Haldane gap Ah which fa- 
vors a disordered phase. For half-integer spin systems, 
the system would be ordered everywhere except exactly 
at the maximally frustrated point where it would be crit- 
ical. 

This behavior predicted by TSDMRG may be found 
analytically by applying the Haldane 0, HJ mapping of 
the Heisenberg model to the non-linear a-model. I write 



f2, 



(3) 



where Qji.Sjl\ttji) = S\Qji) is a spin-coherent state, fiji 
is the local Neel field, r\i is the sublattice modulation, 
m = ^ and 1 describes the ferromagnetic fluctuations 
about the Neel order and Vq is the unit cell volume. Keep- 
ing only up to second order terms in m^;, the Hamilto- 
nian|T|) becomes, 

H' = I dxip^Mx)) 2 + X^m^x) - (Jj_ - 2J d ) x 

n/(a;).n /+ i(a;) + (J± + 2J d )mi(x).m l+1 (x)],(i) 

where p s = -m^^jfVjVjixj - x r ) 2 and x" 1 = 
^jj' 0- ~ VjVj') are respectively the spin stiffness 
and the inverse susceptibility along the chains. As usual, 
one must add the Berry phase term, 



S B = S y^ y r]jUJ [hji] +^2 dxmi(x) 
ji 1 J 



diii (x) 
dt 



x n,(a;).(5) 



Then, writing the total action resulting from H 1 and 
Sb and performing the integration over the fields nii(x) 
with the constraint hi(x).mi(x) — yields, at the maxi- 
mally frustrated point, the following effective action: 
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FIG. 1: TOP: Gaps for S = 1 (left) and S = § (right) for 
J± = 0.2 and for Jd ranging from Jd — 0.06 (bottom) to Jd = 
0.1 (top) with the step 5J d = 0.02. BOTTOM: Extrapolated 
gap for S — 1 and for the same values of Jd as above. 



where K^f = Swx^ 1 + (J± + 2Jd)fiw+i- Hence, to the 
leading order, at Jd = J±/2, the chains are only coupled 
by terms originating from the Berry phase on each chain. 
Therefore as found numerically [J], unlike the magnetic 
case where 5 does not play any role, the physics of the 2D 
system near the maximally frustrated point will strongly 
depends on the ID physics, i.e., whether 5 is integer or 
half-integer as found numerically. 

For relatively small J±, the eventual magnetic order 
parameter is roughly m cx \fj± for Jd = ll| and m 
will get even smaller when Jd ^ 0. For this reason, the 
critical behavior of S = \ systems is very difficult to 



study. I therefore choose 5 — 1 and S = 4 for which m 
is larger and could be extrapolated to from lattice sizes 
not too large. The simulations were done onLx (L + 1) 
systems with L = 8, 12, 16, 20 and 24. I use the periodic 
boundary conditions (PBC) in the direction parallel to 
the chains and open boundary conditions (OBC) in the 
transverse direction. I kept a maximum of ni\ — 243 
states during the first step. I targeted the spin sectors 
S z = 0, ±1,±2. The maximum truncation error during 
this step was p\ — 9 x 10~ 6 . During the second step I 



kept a maximum of 777,2 = 100 states. During the second 
step I targeted two states S z = 0, 1. The truncation error 
was about 5 x 10~ 4 in the magnetic phase. It dropped 
to 6 x 10~ 7 in the disordered phase. The TSDMRG is 
thus at its best in the vicinity of the transition and in 
the disordered phase. 

In Fig[H I show the finite size gap A for 5 — 1 and 
5 = §. In all cases, J± = 0.2 and Jd is varied from 0.06 
to 0.1. The 5=1 and 5 = | show a striking difference. 
For 5 = 1, starting from 0.06 where the system is gap- 
less and decay to faster than a gap opens around 
Jd = 0.073. This gap opening signals a transition from an 
AFM to a magnetically disordered phase. On the other 
hand, the 5 — § system remains gapless for all values of 
Jd- The fact that 5 = | remains gapless makes the tran- 
sition at the maximally frustrated point quite difficult to 
study because long-range order vanishes only at the criti- 
cal point an beyond this point, the Neel phase (ir, 0) sets 
in. This is complicated by the fact that for finite sys- 
tems with OBC, the maximally frustrated point is not 
exactly at Jd = J_i_/2 Hence, the critical behavior 
could only be analyzed for 5=1. I show the in FigQ] 
extrapolated gap as function of Jd for 5 = 1. A vanishes 
at Jd c ~ 0.073. Taking as granted the large TV predic- 
tion that the transition is of second order, I extracted 
the critical exponent zv — 1.205 (A cx (Jd — Jd a ) zv )- 
Hence, Hamiltonian {1} does not belong to the universal- 
ity class of the classical 0(3) Heisenberg model for which 
zv = 0.7048 ±0.0030. 

The center-to-end spin-spin correlation function Cl = 
^(Sl/ 2 ,l/2+iSl : l/2+i) is shown in FigH C L is also de- 
pendent on 5. For 5 = 1, Cl first decays slower than 
L _1 and extrapolates to a finite value for small Jd, then 
decays faster than L^ 1 for larger Jd, where it extrap- 
olates to 0. For 5 = | for all values of the coupling 
studied, Cl extrapolates to a finite value. The behav- 
ior of the magnetization m = v§C^ shown in Fig[2] is 
consistent with that of A. For 5=1, starting from the 
AFM phase for Jd = 0.06, m vanishes around Jd ~ 0.076. 
The best fit to data yields the exponent m cx ( Jd a — JdY ' , 
(3 = 0.3653 which is in good agreement with that of the 
classical Heisenberg model /3 = 0.3639 ± 0.0035. By con- 
trast, to for 5=| extrapolates to a finite value for all 
values of Jd ^5 0.1. This result, with that seen for A, 
contradicts the large N prediction that there is a VBC 
phase in this regime for 5 = | as well. Ultimately, I 
find that long-range order vanishes when Jd ~ 0.11 for 
all L studied and, immediately, the systems jumps to the 
Neel phase with Q = (ir,0). The critical behavior for 
5=| was quite difficult to study. This is because, close 
enough to the critical point, for a fixed Jd, starting from 
the Neel phase with Q = (tt,tt), the system evolves to 
the Q = (7r,0) Neel phase at larger L. For this reason, 
the extrapolations cannot be reliably made. 

A previous Monte Carlo simulation 0] for a Heisen- 
berg model with 5 = 1/2 on the CetV^Og lattice with 
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FIG. 2: TOP: End-to-center correlation C L for S = 1 (left) 
and 5=1 (right) for Jx = 0.2 and Jd ranging from Jd = 
0.06 (bottom) to J d = 0.1 (top) with the step (5J d = 0.02. 
BOTTOM: Extrapolated magnetization for S = 1 and for 
the same values of Jd as above. 



found v — 0.685 and (3 = 0.345. These exponents are 
in agreement with the classical Heisenberg model within 
statistical errors. It is not definitive whether the QMC 
results are in contradiction with my results. In the QMC 
study, following the assumption that critical exponents 
are the same in either side of the transition, v was cal- 
culated from the ordered side through the spin stiffness 
Ps fx (9c — g) zlJ ■ One of the standard results of the clas- 
sical critical phenomena is that the critical exponents 
of the region above the transition are identical to those 
below the transition. It is not however obvious that 



this result extends to quantum phase transitions. From 
the QMC results and mine, it seems that the exponents 
found for model ^ violate this LGW classical behav- 
ior. The exponents in the AFM phase verify the relation 
(d + z — 2 + rj)v = 2(3. Since d = 2, if I assume that the 
critical exponent on either side of the transition are iden- 
tical, I may use zv found in the disordered case in this 
relation. Then, it is clear the above relation between crit- 
ical exponents is not satisfied given that r\ is predicted to 
be small, r\ — 0.0033 for the classical Heisenberg model. 
This violation is possibly related to a dimensional reduc- 
tion at the critical point. At the maximally frustrated 
point the bond strength in the transverse direction van- 
ishes. Hence, the system is effectively ID. It could be 
that this ID physics emerges at the critical point. An 
alternative possibility which is consistent with the small 
difference find in J& c from A and from to, is that the 
transition is of first order and there is a narrow region of 
coexistence of the two phases. 
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No. DMR-0426775. 
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